The Schrodinger Equation for Open Systems 
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An universal exact description of kinetics of open quantum systems in terms of random wave 
functions and stochastic Schrodinger equation is suggested. It is shown that evolution of random 
quantum states of an open system is unitary on average, and this implies validity of the optical 
theorem for any inelastic scattering. 



PACS numbers: 05.30.-d, 05.40.-a, 05.60.Gg 



0^ 
O 
O 



Oh; 

a ■ 
^3 : 



o 

rn 
O 

ON 

o 



- 

X 



1. As it is known, one not always can attribute to a 
quantum system some definite wave function. Instead a 
density matrix is necessary [ij. Anyway if the matter 
concerns an open system S which is a part of more com- 
pHcted closed system S+W (where W is, for instance, 
"thermal bath" or "the rest of world"). In such a case 
even density matrix of S is not quite definite because 
it is impossible to write out an exact general equation 
for it. In practice one is supposed to manage with "ki- 
netic equations" deduced either from the von Neumann 
equation for the whole system S+W or from specula- 
tive reasonings 0, H, 0, Sometimes irreversible terms 
there are supplemented by noise terms to be averaged 
later. Somehow or other, usually one at first constructs 
approximate equations and then constructs a method of 
their approximate solving. 

Meanwhile one may avoid such duplication of approxi- 
mations if recruit "stochastic representation of quantum 
interactions" (see also The latter produces a 

closed stochastic equation for S which is extremely sim- 
ple, intuitively obvious and at the same time formally ex- 
act: if random sources in it reflect, by definite universal 
rules, internal dynamics of W , then such sources act ex- 
actly as actual dynamic interaction of S with W . Impor- 
tantly, these random sources introduce at once noise and 
irreversibility as a statistical effect of the noise. Therefore 
approximation is required only once, under calculation of 
statistical averages. Related mathematical problems are 
similar to problems of theory of oscillations and waves in 
linear media with fiuctuating parameters [ld| . 

A purpose of the present article is to emphasize pos- 
sibility to reformulate the approach suggested in Q (see 
also [3, [Tl|, [1^) in terms of stochastic Schrodinger 
equation and thus generalize to open systems the con- 
cepts of pure quantum state and wave function. Corre- 
spondingly, one can apply to open systems, in an elemen- 
tary statistical sense, the concept of unitarity of evolution 
which naturally coexists with irreversible evolution of av- 
eraged wave function and density matrix of S . As the 
result, the problem of calculation of its averaged (that 
is true) density matrix can be substantially simplified. 
At that, an important roles is played by average wave 
function. 

2. We will start with brief summary of the stochastic 



representation, following ^] (with not great difference of 
designations) and Let us divide the whole system 

Hamiltonian into three parts and write a part Hint re- 
sponsible for the interaction as a bilinear form: 



H = Hs 



Hv 



Hi, 



H„u = E, W, , (1) 



where operators Sj Wj act in different Hilbert spaces 
of S and W , respectively. Such decomposition of Hint 
frequently arises from very physical nature of the interac- 
tion and in any case is feasible in formal sense Then 
for any operator O introduce Liouville super-operator 
C{0) and Jordan super-operator 11(0) as follow: 

C{0)A = ^{AO- OA) , U{0)A = i {AO + OA) , 

where A is arbitrary operator. Thus C{H) is Liouville 
super-operator in the von Neumann equation p = C{H)p 
for density matrix p of the whole system S+W . One 
can easy verify that correspondingly to ([1]) it expands to 

C{H) = C{Hs) + C{Hw)+ (2) 

+ Y,['^{s,)\i{w,) + n{s,)c{w,)] 

Therefore, density matrix ps{t) — T^^w p{t) of S can 
be expressed by 

ps{t) ^lrwm[So^s{t')dt''\p{{)) , (3) 

where fep is chronologically ordered exponential, 

Cs{t) ^ C{Hs) + [x,{t) C{S,) + y,{t) n(5,) ] , (4) 

and Xjit) and yj{t) are time- varying super-operators. 



Xj{t) 



-C(Hw)t 
-C(Hw)t 



Il{W,)e 



C(Hw)t 
,C(Hw)t 



whose action is completely localized in the W 's Hilbert 
space. The latter circumstance, if considered from the 
viewpoint of S , allows us to treat xj (t) and yj (t) in 
^ and ([3]) merely like arbitrary commutative (scalar) 
time-dependent variables. In other words, like "noises" 
or random processes. At that, the trace operation in ([3|) 
plays the role of statistical averaging in respect to the 
random processes Xj (t) and yj (t) . 
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In order to define efi^ective statistical properties of 
Xj (t) and Uj (t) in a simple unambiguous way, it is con- 
venient to assume that before some initial time moment 
to , e.g. ^0 = ; the subsystems S and W were non- 
interacting and statistically independent one on another: 
p{to) = Pg^Vw''' • Then one can write 



ps{t) = {Rit)) , m = Csit)R{t) , 



(5) 



where R{t) is "random density matrix" of S , which 
obeys "stochastic Liouville equation" with initial condi- 
tion -R(to) = Ps""* I angle brackets symbolize the 
averaging: 

(...) = Trvy ... Pw' 

On right-hand side here Xj{t) and yj(t) are treated as 
the above super-operators while on left side as effective 
random processes. This is just the statistical definition of 
these processes. Evidently, it means that characteristic 
functional of random processes Xj (t) and yj (t) can be 
presented by expressions 

( exp y ^ [ u, (t) X, it) + f, it) y, (t) ]dt} ^ (6) 
Cw{t) dt 



Trw exp 



(in) 

Pw 



where Uj{t) and fj{t) are arbitrary probe functions, 

Cwit) = C{Hw) + E[^j(t)^iWj) + f,it) C{W,)] , 

the integration begins at t = to and by request one may 
move ^0 away to — oo . 

According to ^ ©, the processes x{t) represent per- 
turbation of 5 by W as directly prescribed by their 
Hamiltonian interaction. Analogously, y(t) represent 
back action of S onto W . In view of these observations, 
we can expect that x{t) are quite usual random processes 
while, in opposite, y(t) are rather unusual ones. Indeed, 
at u{t) = the super-operator Cw{t) reduces to a Li- 
ouville super-operator, and then from ^ it follows that 
{y{ti) ... y{tn) ) — ■ By the same reason, as we can see 
from (111), {y{t) x{ti) ... x{Tra)) = if i > maxfc(Tfe) . 
But generally {x{ti) ... x{Tm) y{ti) ... y{tn)) / if 
maxfc(r/j) > maxk{tk) . Hence, y{t) show their worth 
in conjunction with x{t) only. At that, cross correla- 
tions of y{t) with x{t) represent response of W to its 
perturbation by S and eventually introduce to evolu- 
tion of ps{t) dissipation and irreversibility. The condi- 
tion maxfe(Tfc) > maxfc(ifc) comes from the chronolog- 
ical ordering in ^ and expresses the causality princi- 
ple. Details and generalizations of such "stochastic rep- 
resentation of dynamic interactions" Q can be found in 

[i,0,i,ilii|,[i2,lil,lil. 

3. Now let us go to original theme of this paper. 
If initially (before interaction between subsystem un- 
der interest S and "the world" W was switched on) S 



was in a pure state, pg"^ = |\I'(™))(vI/(*") | , then solution 
to equation (O is merely R{t) = |\I'(t))(*(t)| , where 
'i'{t) satisfies "stochastic Schrodinger equation" : 



^-[Hs + y] w,{t)S,]^{t) 



(7) 



dt n 

with initial condition '5 (to) = vl/f*"' and 

Wj{t) = Xj{t) -I- ihyj{t)/2 

being interpreted as complex random processes. In such 
sense, S at all times remains in a "random pure state" . 
Consequently, at arbitrary initial density matrix we can 
make its diagonalization, p^™' = J2a I^q"''} Pa (5'q"''| , 
and then write 



(8) 



where each of ^'q-(<) satisfies ([7]) under initial condition 

^0(^0) = ^a""* , while the weights Pa stay constant. 

Because of complexity of Wj (t) an evolution described 
by equation ([7]) is not unitary: 



dt 



{^a{t)\^0{t)) = ^ 



Nevertheless, it is unitary on average: 

I {{^am.it)))^ (^ y,it){^a{t)\s,\^m) = 

This statement directly follows from the above mentioned 
statistical properties of y{t) (that is factually from the 
causality principle). As the consequence, 

( (^P„(t)|vi/^(o) ) = (*i'")|*{;"^) = Sap , (9) 

and one can say that the "random pure states" all the 
time are mutually orthogonal on average. 

Let us suppose that the interaction by its nature is 
localized at time, that is has character of scattering pro- 
cess. In such a case, random time- varying state of S can 
be written as sum 

\^it)} ^ 1*0(0) + 1*^(0) , 

where ^'^(t) represents free evolution of S while ^'^(t) 
is contribution of the scattering, so that ^''*(io) = 0. 
Then, due to the unitarity on average, as well as unitarity 
of the free evolution, the equality 

{{^m{t))) = {^'m"{t)} 

takes place. It implies that 

2Re(*"|(*")) + ((^'l^-*)) = (10) 

This is generalized "optical theorem" which is valid for 
arbitrary non-elastic scattering. 
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This result demonstrates that the average wave func- 
tion (\E'(i)) contains important information about an 
open system. It should be underlined that ( ^'(i) } is al- 
ways unambiguously defined by equation ([7]) as soon as 
random processes w{t) are defined by (O. Therefore, 
the average wave function is not less legitimate than the 
density matrix ps{t) . 

Notice also that if 'w{t) in ^ are stationary random 
processes then average wave functions {'^^{t)) in (fTO|) 
and = -I- include those energies 

and frequencies only what already are present in '^'^{t) . 
Therefore non-elastic component of the scattering wholly 
comes from fluctuations ^{t) — {"^{t)) , although the lat- 
ter contribute to elastic component too. In view of this 
circumstance, relations ([9]) or (fTO|) imply inequality 

prnel^^^ < I - { {t)) \ (t)) ) , 

with V (t) being total probability of non-elastic scat- 
tering. In case of a scattering center in one-dimensional 
conducting channel one can obtain additional inequality 

J,^nel < [ll|, . 

4. The stationarity condition realizes, for instance, 
when pI^'' cx exp (— iJvt^/T) , that is subsystem W ini- 
tially was thermodynamically equilibriuni (all the more 
when W is a thermostat). Then 0, 0, B El, E, [11 

{w*{T),Wmm = Kjmir) , 
{Wj{T) ,WmiO)) = Kjjni\T\) , 



-f exp(?iw/r)e- 
1 + exp{nuj/T) 



where aj^i^^) is a non-negatively defined real- valued 
spectrum matrix. This information about the random 
sources w(t) is sufficient in the framework of the Born 
approximation or "one-loop" and "ladder" approxima- 
tions and their analogues (such as "Bourret approxima- 
tion", "Kraichnan approximation" , etc. [lo|). Moreover, 
if w{t) possess Gaussian statistics ( W is "Gaussian ther- 
mostat") then the spectrum ajmii-o) contains exhaus- 
tive information about them, and principally one can ex- 
actly find any average value. In all that cases, organi- 
zation of the correlation matrix Kj„i{T) automatically 



ensures |ll|, [ij] balance of noise and dissipation (ra- 
diation and absorption of energy by W ) in agreement 
with the fluctuation-dissipation theorem. If w{t) are 
non-Gaussian then they obey also "generalized non-linear 
fluctuation-dissipation relations" 14| . 



5. Thus, the traditional approach to open sys- 
tems based on model kinetic equations (or equivalent 
"Langevin equations" ) has an alternative in the form of 
exact linear Liouville and Schrodinger stochastic equa- 
tions supplemented by exact universal expression for 
characteristic functional of random sources (fluctuating 
parameters) in these equations. At this alternative ap- 
proach it is possible to exploit concepts of wave functions 
and (pure) quantum states, their superposition and uni- 
tary evolution. Correspondingly, mathematical difficul- 
ties concentrate at statistical averaging of wave functions 
and density matrices over the sources. This allows to use 
many achievements of the theory of linear systems with 



ory c 

randomly varying parameters [lO|. In this respect it is 
useful to notice that non-Gaussian statistical models of 
fiuctuating parameters sometimes are much better suit- 
able for exact calculations than Gaussian ones 
Therefore development of simple Hamiltonian dynamical 
models of non-Gaussian thermostats is of great interest. 
From the other hand, as was shown b y e xample appli- 
cations of the new approach [1, 0, [ll|, 12 1, quite useful 
results can be obtained even without full statistical in- 
formation about the sources. 
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